Electronic Structure and Lattice dynamics of NaFeAs 
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The similarity of the electronic structures of NaFeAs and other Fe pnictides has been demonstrated 
on the basis of first-principle calculations. The global double-degeneracy of electronic bands along 
X-M and R-A direction indicates the instability of Fe pnictides and is explained on the basis of a 
tight-binding model. The de Haas-van Alphen parameters for the Fermi surface (FS) of NaFeAs 
have been calculated. A Qa/ = (1/2, 1/2,0) spin density wave (SDW) instead of a charge density 
wave (CDW) ground state is predicted based on the calculated generalized susceptibility x(q) an d a 
criterion derived from a restricted Hatree-Fock model. The strongest electron-phonon (e-p) coupling 
has been found to involve only As, Na z-direction vibration with linear-response calculations. A 
possible enhancement mechanism for e-p coupling due to correlation is suggested. 
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I. INTRODUCTION 

The recent discovery of superconductivity in 
LaFeAs[Oi_. E F. E ] 1 with the remarkable critical tem- 
perature T c of 26 K has stimulated intense experimental 
and theoretical research activities towards exploring new 
superconducting systems and understanding the relevant 
physics. So far, four systems represented by RFe(As,P)0 
(R=La,Ce,Pr,Nd), (1:1:1:1), (Ba,Sr,Ca)Fe 2 As 2 , (1:2:2), 
(Li,Na)FeAs (1:1:1) and Fe a (Te,Se,S) (1:1) have been 
discovered. They all share the similar structural 
character of two- dimensional Fe-As layers formed by 
face sharing FeAs4 tetrahedra sandwiched by ionic 
layers such as ROi_ x F x , Bai_ x K x , Lii_a;,Na x , or by 
van der Waals gaps in the (1:1:1:1), (1:2:2), (1:1:1), 
(1:1) systems, respectively. The stoichiometric systems 
can be doped with electrons or holes by appropriate 
substitutions in the ionic layers. Without doping most 
of the parent systems show universal spin density wave 
(SDW)-like long range magnetic order and a tetragonal 
to orthorhombic/monoclinic structural phase transition 
at low temperature.— The doping supresses both the 
magnetic and structural phase transitions and leads to 
superconductivity. The (1:1:1) system seems to be an 
exception, in which no phase transition was observed in 
early experimentsi^' 5 -'^ 

Density functional studie a 7 ' 8 ' 9 ' 10 ' 11 have shown that 
the electronic states around the Fermi level are of pre- 
dominant Fe-3d character. The minor splitting of these 
d bands is mainly due to the weak tetrahedral ligand 
field of As. The universal feature of nesting conditions 
for parts of FS is generally thought to be the origin of 
the magnetic and structural phase transitions. The miss- 
ing transitions in the (1:1:1) system call for an expla- 
nation. So far metric details of the FS for all systems 
have not been addressed, which would be important to 
provide quantitative tests to the present understanding 
of Fe-based superconductors when compared with dHvA 
data. Furthermore, the role of phonons remains impor- 
tant even though earlier calculation o 12 ' 13 which show a 
rather weak electron-phonon (e-p) coupling, as the im- 
portance of magnetophonon coupling is revealed in a re- 



cent calculation. 10 In this work, we investigate the elec- 
tronic structure and lattice dynamics of NaFeAs, because 
of the aforementioned behavior and its simple structure. 



II. COMPUTATIONAL DETAILS 



For the numerical calculations the experimental crys- 
tallographic data of NaFeAs at 2.5K 5 - were used. The 
magnetic ground state is not chosen as the starting point, 
because we seek a comparison with the electronic struc- 
tures calculated for other Fe-As based systems. The full- 
potential LMTO metho d 14 ' 15 based on the density func- 
tional theory in the form of LDA+GG A 16 ' 17 is used for 
all calculations. The muffin-tin (MT) radii of 2.249,2.340 
and 3.271 a.u. for Fe, As, Na, respectively, are calcu- 
lated according to the maximum of overlapping Hartree 
potentials of neighboring atoms.— The analysis of states 
is based on these MT radii. A 2-k basis set with tail 
energies of -1.10, -0.1 Ry, respectively, was used for the 
Fe-4s, -4p, -3d, As-4s,-4p,-4d and Na-3s, -3p, -3d states. 
While the Fe-3p, As-3d and Na-2s, 2p states were treated 
as semi-core states, all the others were handed as core 
states. The charge densities and the potentials inside the 
MTs were expanded into harmonics up to lmax=6, while 
those inside the interstitial region were expanded into 
32646 plane waves corresponding to a cutoff energy of 
190.22Ry. A k-mesh of 32x32x28 was used to calculate 
the electronic structure. As the calculation of the gener- 
alized electronic susceptibility x(Q? w ) is sensitive to the 
k-mesh, we have chosen a mesh of 44x44x40, which is 
similar to that used in Refi 7 -. For the linear-response cal- 
culations, a coarse mesh of 12x12x12 and a dense mesh 
of 36x36x36 were used. The latter was used to attain the 
accurate integration weights in the linear-response calcu- 
lations. A q-mcsh of 6x6x4 was adopted to produce 
30 independent phonons. Except for these parameters, 
all of the other parameters are the same as those in the 
electronic structure calculations. 
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FIG. 1: Band structure calculated from FP-LMTO method, 
the black dot indicating a representative state, ipz,i2, for the 
steep band. The curves in the right part of the figure are solid 
for total DOS, blue broken line for Fe-3d, red broken line for 
As-4p, and black broken line for Na. 



III. RESULTS AND DISCUSSION 
A. Electronic structure: General features 

The calculated FP-LMTO band structure and the elec- 
tronic DOS for NaFeAs are shown in Fig. Q] It is clear 
from the projected DOS curves that within the energy 
window of Ef ± 2 eV, the electronic states are of pre- 
dominant Fe-3d character which determines the physi- 
cal properties of the system. By using the TB-LMTO 
method^, we have calculated the hopping integrals de- 
fined as 



(i) 



where \i) indicates a orbital at site i and H is the LDA 
Hamiltonian. The three most important hopping in- 
tegrals for adjacent Fe ions (2.791A) are H d d =- 
0.284, H dx2 _ y2 _ dx2 _ y2 =0.230, and H dxy ^ 2 _=-Q.l^ 
eV, while the most significant hopping integral for the 
next nearest Fe neighbors (3.947A) is H d ^ 2 _ 2 - d ^ 2 _ 2 =- 
0.036 eV. Below a small energy gap the hybridiza- 
tion region of covalent Fe-As and partially covalent 
Na-As bonding sets in. The calculated dominant 
hopping integrals between As and Fe (2.428A) are 



H Py - dyz =0.816, H p ^_ 



=-0.620, H n 



=-0.428, 



H Py - d;j 2 i =0.427 eV. Though the hopping integral be- 
tween Na-s and As-s,p orbitals are quite large (e.g. 
# s _ s =1.256 eV, ff s _ s =1.218 eV for the Na-As distances 
2.984 and 3.107A, respectively), the rather small pro- 
jected DOS of Na in the energy range of interest leads to 
a much weaker Na-As covalent bonding as compared to 
Fe-As and Fe-Fe. The strong Fe-As covalency determines 
the bonding in the 2D Fe-As layer, while the partial ionic 
bonding between Na-As determines the packing of these 
layers. It in important to emphasize that the covalency 
between Na and As is neglected in many studies, due to 
the choice of a smaller MT radius for Na, which assigns 
some electronic states to interstitial electronic states. 
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FIG. 2: The left figure shows the structure of the Fe-As layer 
in NaFeAs with ti, i=l,2 representing the primitive lattice 
vectors of the crystallographic unit cell; t[ are the primitive 
vectors implied by the electronic band structure. The right 
figure shows the corresponding reciprocal lattice vectors and 
BZs. 



The electronic structure as shown in Fig. [T] is very 
similar to that of LiFeAs calculated by other methods i^iSii 
Though the detailed comparison of electronic structures 
of different systems is difficult due to limited data and dif- 
ferent choices of symmetry directions, the overall agree- 
ment of the energy spectrum around the Fermi level is 
obvious. 



B. Electronic Structure: Double Degeneracy 

An interesting feature which has not received much 
attention is the double degeneracy of all bands along X- 
M and R-A directions. Nekrasov et al^- have mentioned 
that there are two degenerate bands along the X-M direc- 
tion in the (1:1:1:1) system, however, without any further 
discussion. As the symmetry of the X-M line is identical 
to that of the R-A line, we here concentrate on the X-M 
line. The symmetry groups at X and M points are D2h 
and -D4/1, respectively. The D2/1 symmetry at X and the 
Civ symmetry of points along X-M do not enforce any 
degeneracy of the bands, since Din and Ci v have only 
ID representations. The D±h symmetry at M, however, 
results in double degeneracy for bands belonging to the 
Eg or E M representations. From M to X, this degeneracy 
will generally be reduced due to symmetry lowering, e.g. 
Eg — > A2 + B2 ■ Obviously, the global double degeneracy 
of bands along the X-M line is not caused by the known 
symmetry. Considering the global property of this de- 
generacy, a hidden symmetry needs to be investigated. 

The predominant Fe-3c? character of all bands within 
Ef±2 eV indicates that these bands originate from pure 
Fe lattice. It is obvious from Fig. [2] that each of the two 
As square lattices projects onto the center of the other 
lattice, each of which is crystallographically equivalent to 
the Fe lattice. Therefore a tight-binding energy disper- 
sion for a Fe-d band can be derived from a 2 x 2 secular 
matrix as follows, 



Hn — £ H12 
H21 H22 — £ 



= 0, 



3 



where the indices 1,2 indicate Fel and Fe2 atoms in the 
unit cell, respectively. The on-site energy Hu = i? 22 of 
the chosen orbital is calculated up to the second nearest 
neighbors, and is given for the X-M line, by 

H n =a d + (3f (2 cos(27rx) - 2) - 4/3^ cos(2ttx), (2) 

where a^jPf,^ are the on-site energy, nearest neigh- 
bor, and second nearest neighbor hopping integrals be- 
tween Fel-d orbitals, respectively. The parameter x € 
[0,1/2] specifics a k point along the X-M line with 
k = xgi + l/2g 2 , where gi and g2 are the reciprocal 
lattice vectors corresponding to ti and t 2 , respectively. 
An interesting result is that 7J 12 = H 21 = for k points 
on the X-M line. In the tight-binding approximation the 
conditions of Hu = H22 and Hyi — = result in a 
doubly degenerate band dispersion along the X-M line, 
i.e. E(k) = Hu. 

In fact, if only Fe bands are concerned, a smaller 
unit cell as shown in Fig. [2] for the Fe sub-lattice, with 
the unit cell vectors, t l = l/2(ti — 1 2 ), t 2 = l/2(ti +t 2 ), 
t 3 = t3 can be found, which results in a double sized 
first BZ with the X-M line as an internal line. When 
the reduced zone scheme is used to represent the bands, 
the bands in the large BZ will be folded into the smaller 
BZ as shown in Fig. [2] by the reciprocal lattice vector. 
Accordingly the degeneracy of the boundary states will 
be doubled, as can be easily understood from the one- 
dimensional superstructure case i 22 i 23 

It is obvious that the band-folding explanation can 
not be directly applied to the As bands or to the Fe, 
As hybridized bands, since the As sublattice can not be 
described by t { ,i = 1,2,3. The two As atoms in the 
unit cell are connected by a glide plane operation. To 
understand the double degeneracy of these bands along 
X-M, we analyze the Hamiltonian matrix constructed by 
tight-binding approximation. We do not attempt to solve 
the whole secular problem which is a 28 x 28 matrix, if 
all Na-3s, As-4s, 4p, and Fe-4s, 4p, 3d orbitals are consid- 
ered. Instead we investigate the most representative and 
important subblock by exploiting the symmetry. Guided 
by the electronic structure shown in Fig. [TJ we consider 
in the following the As-p and Fe-d orbitals only. By ac- 
counting for the C 2 „ symmetry of the X-M line, we choose 
the £?2 subblock as an example. As only p z and d xz be- 
long to the £?2 irreducible representation, we attain a 4 x 4 
matrix secular problem: 



/Hu 


-e 










H n 


- £ H h 








H32 




-e 


\H 41 








H 33 — e 



where the Hamiltonian matrix elements in both the rows 
and columns are arranged in the order of Fel, Fe2, Asl 
and As2 (Fig. ^ , respectively. Obviously the upper 2x2 
diagonal subblock corresponds to the Fe sublattice, while 
the lower diagonal subblock corresponds to the As sub- 
lattice. Though the As sublattice has a different symme- 
try as compared to the Fe sublattice, their tight-binding 



Hamiltonian matrix have the same diagonal structure. 
Our calculations based on the structure as shown in Fig. 
[5] indicate that H 43 = H^ 4 = 0, while H 33 = H44 cal- 
culated up to the 2nd nearest neighbors of As is given 
by 

H 33 =a p +/3{ (2cos(27nzi) - 2) - 4/3f cos(27re), (3) 

where a p , 0^,02 are the onsite energy, nearest neigh- 
bor, second nearest neighbor hopping integrals between 
Asl-p z orbitals, respectively. The similarity between the 
Fe, As subblocks of the secular matrix suggests that the 
bands of the As sublattice should also be doubly degen- 
erate along the X-M line. When the Fe-As interactions 
are taken into account, the situation becomes more com- 
plicated, because there are non-zero off-diagonal matrix 
elements. Under the aforementioned approximations, the 
matrix elements H 3 i — i?* 3 = and H42 = H 2A = are 
obtained. 7J 32 is found to be equal to the complex con- 
jugate of i?4i, i.e. H32—H41, and is given by 

H 32 = P P d (1 + e~ i2 ™) + fi\ d (-2 - 2e-' i2 ™) , (4) 

where (3 p d, /3f d , are the nearest neighbor, second nearest 
neighbor hopping integrals between Asl-p z and Fe2-d xz 
orbital , respectively. Simple algebraic operations of the 
4x4 secular matrix equation lead to the following reduced 
problem: 

((Hu ~ e)(H 33 ~ e) - \H 32 \ 2 ) 
x ((Hu - e)(H 33 - e) - \H 41 \ 2 ) = 0. (5) 

Since ff 32 = H\ x , Eq. [5] which is of fourth order is re- 
duced to two identical equations, resulting in two identi- 
cal sets of roots, namely, 



e= 11 2 33 ±-^(Hu-H 33 f+A\H 3 2\*. (6) 

This result provides an explanation for the double de- 
generacy of bands along X-M line. As can be seen from 
Eq. [6] if the coupling between Fe and As sublattices is 
removed, i.e. -ff 32 = 0, the situations of pure Fe and As 
lattices are recovered. Within our tight-binding model, 
the occurrence of zero matrix elements and their posi- 
tions are essential to the double degeneracy. Our cal- 
culations show that this is a combined result for the site 
symmetry of Fe, As and the special direction X-M, which 
makes the terms e ini2nx e m2 * H i:j , e ini27TX e in '^ H i:j with 
\n 2 — n' 2 \ = odd, always occur in pairs and thus cancels 
with each other, where niti +n 2 t 2 indicates the position 
of the target atom for the hopping integral of interest. 
This conclusion holds even when the interlayer hopping 
terms are included, because t 3 is perpendicular to k along 
the X-M and R-A directions, and the algebraic structure 
of the secular matrix remains unchanged. Therefore, for 
the bands along X-M or R-A the system can be described 
with the smaller unit cell as shown in Fig. O Granted 
other conditions remain unchanged, the reduction of a 
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FIG. 3: The calculated FS from FP-LMTO method are shown 
on the top, the first three objects of hole nature center at F 
with F-Z as their center axis. The last two objects of elec- 
tronic nature center at M with M-A as their center axis. The 
arrows on the surface show both the magnitude and direction 
of the velocity field mapped onto the FS. The first BZ used 
throughout this work is shown below the FS. 



unit cell implies a higher symmetry in the sense that 
more translational operations are allowed. The global 
degeneracy of bands along X-M and R-A directions sim- 
ply infers an instability of the system. 



C. Electronic Structure: Properties of Fermi 
Surface 

The calculated FS shown in Fig. [3J is composed of 
a small ovoid ball, two large tubes with hole character 
around T and the other two tubes with electron charac- 
ter around M, respectively. Qualitatively these features 
agree quite well with those of other systemsiZ>2»2S A sim- 
ple visual inspection reveals possible Fermi surface nest- 
ing between the hole tubes around T and electron tubes 
around M corresponding to a vector of (1/2,1/2,0). The 
calculated group velocities of electrons/holes are mapped 
onto the Fermi surface as shown in Fig. [3J where the 
group velocity is defined by, 



v(kj) = ~V£(kj). 



(7) 



The calculated Fermi velocities of electrons are much 
larger than those of the holes (see Fig. [3]). The direc- 
tions of the velocities suggest that the movement of most 
electrons and holes are confined to the Fe-As layers. Only 
a few holes which are located at both ends of the ovoid 
ball move perpendicular to the layers. This can also be 
seen in Fig. [TJ where a steep band crosses the Fermi level 
along the T-Z direction. The property of this band will 
be discussed later. So far, there have been no experimen- 
tal result to prove the existence and quantitative details 
of the FS for NaFeAs and also for other systems reli- 
able experimental data have not been established. Cal- 



TABLE I: Calculated Metric Details of Fermi Surface and 
Band Masses 



Orbit" 



F eaie (Tesla) 



rh 
- V 

rh 
- z 

Zh 
■ , z 

rh 
.- z 

Zh 
■ - z 

15 Ml 
16-At 



627.3 
1452.2 
1692.9 
1666.0 
2078.6 
2155.4 
1836.7 
2820.9 
1392.2 
2116.3 



-.473 
-.776 
-1.940 
-1.954 
-1.480 
-1.415 
0.918 
1.673 
1.108 
0.810 



"Orbit is symbolized with band index k norma i > where k indi- 
cates the center of the orbit; F ca [ c indicates the calculated area of 
an orbit with specified center and normal;Band masses are given 
in unit of free-electron mass. 



culations emphasizing the correlation effect^ predict a 
large scattering rate and short life time for states around 
the Fermi level, and thus disappearance of hole pock- 
ets at r even in the case of a doped (1:1:1:1) system. 
The issue here is whether LDA has considered in a suf- 
ficient way the correlation effects. For a further test of 
different approaches, we have also calculated the met- 
ric details of the Fermi surface as summarized in Table 
HI The band masses shown in Table H] are calculated as 
mcaic = dA/dEi 25 i 26 It needs to be pointed out that 
the calculated Fermi surface may split due to a SDW of 
the actual ground state. However, as this split is not 
too large, the calculated area of various orbits may still 
be valid to be compared with the de Hass-van Alphen 
experiments where available. 



D. Electronic Structure: SDW or CDW instability 

The above discussions have already shown that the 
electronic structure of NaFeAs shares many universal fea- 
tures with the Fe-based superconductors. To further in- 
vestigate the instability of the electronic system of the 
paramagnetic state, we have also calculated the general- 
ized susceptibility x(q, = 0) which is defined as follows, 



x(q,«) = y fii:kj: /;Ak -' ! 



i0+ 



(8) 



The integrand of Eq. [8] is a integral kernel for construct- 
ing a T matrix which is directly connected to the dielec- 
tric constant including local field corrections >21 

To investigate the singularity of x(qj0)j lts values 
evaluated within the aforementioned dense q-mesh are 
mapped onto different planes perpendicular to the z axis. 
Our investigations indicate that the peak of i?ex(q, 0) 
occurs for (1/2, 1/2, z), but becomes increasingly weak 



5 




t 



' T N 



t 



b 



FIG. 4: Calculated J?e x(Q; 0) mapped onto the plane through 
F with z axis as its normal. The calculated peak value at M 
is 1.14 x 10 11 with the minimum scaled to zero. The right 
figure shows the SDW derived from Q M = (1/2,1/2,0), by 
assuming a linear polarization. 



from z=0 to the boundary, 2=1/2. This finding together 
with the results for /mx(q, 0) agree quite well with the 
earlier calculations for the (1:1:1:1) system!. Fig. [4^ 
shows the function values of Re x(q, 0) on the T plane. It 
is obvious that a peak appears at q, (1/2,1/2,0), which 
implies that the condition for a SDW, V q > 1/xq can 
easily be satisfied at low temperature ) 28 ! 29 with V q be- 
ing the average exchange matrix element. By assuming 
a linear SDW with p(R) = pxcos(q • R), one attains 
a spin structure as shown in Fig. [4}d, where p has the 
dimension of spin per unit volume, R is the lattice po- 
sition for Fe. The SDW as shown in Fig. 0]d breaks 
the original symmetry and supercell of y/2 x V2 sets in. 
This type of SDW state has been shown by Yildirim 11 
with DFT calculations to have the lowest energy among 
the investigated states for LaFeAsO. It is obvious that 
the nearest-neighbor spin coupling along the two diag- 
onal directions is ferromagnetic and antiferromagnetic, 
respectively. This asymmetry may favor the tetragonal- 
orthorhombic transition as observed in other systems. It 
should to be pointed out that the spin dynamics need to 
be considered in a more strict way since the d-electrons 
around the Fermi level are itinerant. 

The peak of Re x(q, 0) only implies the instability 
of an electronic system, which may result in a SDW, 
or other new orders, in particular a CDW with lat- 
tice distortion as discussed by Chan et. al. based on 
a Hartree-Fock approximation (HFA) with screened ex- 
change interaction^ The condition for the occurrence of 
a CDW within the restricted HFA can be written as 



huj 



(2C/ q - V q ) > 



1 

Xq' 



(9) 



where [/ q and r? q are the average Coulomb interaction 
and e-p interaction matrix elements, respectively. If the 
values of £/ q and V q are estimated to be of the order 
of ~ 4 eV, and ~ 0.7 eV as obtained for LaFeAsO 2 ^, 
then a strong e-p interaction and a soft phonon mode for 
the nesting q vector are prerequisites for the CDW to 
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FIG. 5: Phonon dispersion calculated from linear-response 
theory implemented in the FP-LMTO method. 



take place. Earlier calculations*^ 3 - using other methods 
have already revealed a rather weak average e-p coupling, 
however, no discussions have been made in this respect. 



E. Lattice dynamics and e-p coupling 

The calculated phonon spectrum is shown in Fig. [5l 
where overall agreement with results using ultrasoft pseu- 
dopotential method^ 3 - is attained. The main differences 
occur along T-Z, where obviously the authors of Ref. 13 
have chosen a different coordinate for the Z point. In our 
calculations, we have chosen a Z point which lies at the 
boundary of the standard BZ as shown in Fig. [3l which 
agrees also with that of Ref.—. The detailed compari- 
son for the phonon frequencies at V is given in Table |TT] 
by transforming the units from THz in Fig. [5] to cmT x . 
The overall agreement between two different calculations 
is rather good. The largest differences are observed for 
the two highest IR-active modes, which are only 22.15 
and 21.99 cm -1 for the Ai u and E u modes, respectively. 
These values can easily be tested experimentally. As 
shown in Fig. O the phonon spectrum at M does not 
show any obvious softening. The evaluation of the first 
term in Eq. [9]needs more effort. In the original derivation 
of Eq. [9] unit cells with one more atom were not consid- 
ered. All e-p matrix elements ^(k, k + q, q) were treated 
as an average constant 77 . As a rough estimation, if the 
actual e-p matrix elements, ?y(kj, k + qj', qv) can be re- 
placed by an average, then the condition given in Eq. [5] 
can be used to judge the occurrence of a CDW related to 
a phonon qv. The e-p matrix elements ?y(kj, k + qj', qv) 
can be calculated by linear-responce theory as in Ref. 15 , 
where it is symbolized as fl£+qj' kj -- In this work, we use 
the following Fermi surface average to estimate r/ 2 . 



v 2 



N(0) 



EiC^ k /^)^w) (10) 



where N(0) is the density of states at the Fermi level. 
By using Eq. [TUl the first term in Eq. [5] can easily be 
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TABLE II: The calculated phonon frequencies in cm 1 at Y for NaFeAs, the Raman and infrared-(IR) modes are indicated 
with their irreducible representations in D^h- For comparison the values from Ref.— are listed in the last two rows. 

Raman 125.77-£ 9 173.98-Ai g 198.01-Ai fl 202.65-£ 9 227.35-B i9 

IR 181.56-.Eu \9,7.m-A 2 u 254.99-£' li 275.15-A 2u 

Raman llO-K, 176-Ai g 199-Ai 9 187-£ 9 218-Bi g 

IR 170-Eu 183-.42U 233-_B u 253-A 2ll 




a b 

FIG. 6: q dependence of e-p coupling constant, A(q), in the 
first BZ of phonons revealing a peak-like structure; The right 
figure shows the most important Ai g phonon in e-p coupling. 

calculated using the phonon spectrum displayed in Fig. 
[5J We have calculated j-^- for all other modes except the 
three trivial acoustic phonon modes. At q, (0.5, 0.5, 0), 
the largest value of j— — of the 18 phonon branches is 
0.542 eV with v — 3, which is too small to satisfy Eq. 
[5J This value is smaller than that of Cr as estimated 
from the experimental data^ In fact, even when we take 
the sum of all l^k+qj' kj I 2 vaiues & t the Fermi surface 
without averaging, the largest value is 2.08 eV with v = 
3, which is still too small to satisfy Eq. [51 We thus 
exclude the possibility of a CDW with lattice distortion 
related to the nesting vector. The largest two values on 

our q mesh for |^ are 2.99 and 2.73 eV at Z, (0,0,1/2) 
and q, (0,0,1/4), respectively, with v — 6, which are 
still too small to satisfy the criterion of Eq. [9] This 
fact also implies that the largest e-p interaction occurs 
in the interval between T and Z for phonon branch 6 as 
indicated in Fig. [5] To demonstrate this result more 
clearly, we have calculated the phonon mode dependent 
e-p coupling constant defined as follows, 

\ - 1 

where 7 q „ is the phonon line-width as defined in RefJ^. 
The values A q „ for the above two phonon modes are 1.495 
and 1.365, respectively, which are also the largest val- 
ues among all the phonon modes we studied. This con- 
sistence indicates that under the approximation of Eq. 
[TU1 the first term in Eq. [5J is essentially another mea- 
sure of the e-p coupling. To reveal the q dependence 
of A, we have further summed out the phonon branch 



index v of A q „, see Fig. [f^i. As shown in this figure, 
the characteristic peak-like structure of e-p coupling in 
superconductor a 31 ! 32 ' 33 exists also for NaFeAs. However, 
in the case of NaFeAs, the maximum value is 3.19 at 
q, (0,0,1/4), compared to 12.3 in HgA 25.3 in MgB^ 
and 25 in Lii-^BC 32 . As also shown in Fig. [5^, the 
most effective e-p coupling occurs along T-Z. Except for 
the two largest values, the other effective ones are all 
at least 4 times smaller, while all remaining ones van- 
ish. This result is different from those calculated for 
LaFeAsOi^ajF^, where the e-p coupling distributes rela- 
tively evenly among the q points and phonon branches. 12 
In particular, the most important e-p coupling is at T, 
where we found no coupling at all. The reason for this dif- 
ference needs to be clarified in future work since there are 
many similarities between the two systems. An unusual 
result is that the dominating e-p coupling stems from 
phonons involving only As and Na vibrations perpendic- 
ular to the Fe-As layer. In Fig. [6Jd., the A\ g phonon at 
Z which exhibits the largest X qv value is shown. The less 
effective phonons have \ qv values below 0.3, involving 
the Fe vibrations in the X-Y plane. However, even for 
these phonons the dominating component is still based 
on As, Na vibrations along the z direction. This result 
indicates that only very few electronic states around the 
Fermi level are effectively coupled to the phonons. As can 
be seen from the projected DOSs shown in Fig. [T] the 
weights of As and Na states at the Fermi level are rather 
small. By calculating the compositions of the Fermi 
states, we have found that the As and Na show significant 
presence only in the steep band crossing the T — Z line as 
shown in Fig. [T] A representative state ipz,i2 as indicated 
in Fig. Q] is plotted in Fig. [7] It is understandable that 
the d x 2_ y 2 orbital of Fe is less affected by the vibrations 
of As and Na due to the small overlap between the rele- 
vant orbitals. The large lobe of the Na-s,p hybrid is due 
to its diffuse character, which is often treated as intersti- 
tial electrons in other methods. The weak e-p coupling 
for specific phonons as discussed above and in particular 
the rather small percentage of effectively coupled elec- 
tronic and phononic states in the phase space result in a 
small total A=0.194 compared to 0.27 from the pseudo- 
potential calculations^ 3 - and 0.21 for LaFeAsOi-^FaJ^. 
With this value and the calculated logarithmic average 
phonon frequency, o->; og =215.8K, the calculated supercon- 
ducting transition temperature T c based on the Allen and 
Dynes formula T c is close to zero. Even if the multiband 
effect is considered, the experimental T c of 9K is difficult 
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FIG. 7: A representative state for the steep band which cou- 
ples effectively with the Ai g phonon. The largest lobes indi- 
cate the Na-s,p hybrid, the middle sized lobes indicate the 
As-s,p hybrid, the most contracted lobes indicate the hybrid 
of mainly Fe-d x 2_ 2 character. 



prohibiting factor for a CDW, however, the correlation 
effect needs further discussion. It has already been ar- 
gued that the Fe pnictides are also in proximity to Mott 
insulators i 24 ! 37 By using the calculated potential param- 
eter, we have estimated the band width W3d=3.5 eV for 
the Fe-3ci orbital as compared to the Coulomb param- 
eter U ~ 4eV. This implies that NaFeAs should be a 
Mott insulator or more accurately a bad metal consider- 
ing the degeneracy of d bands. 24 When correlation effects 
are considered at the DMFT level 2 ^, very large spectral 
weight transfers from B\ g —d x 2_ y 2 to the Ai g ~d 3z 2_ 1 and 
a small increase for the As-4p orbital occur. Considering 
the established similarity between the two systems, we 
would expect an analogous transfer to take place, which 
means the d x 2_ y 2 component in Fig. [7]will be replaced by 
d 3z 2_i and thus increase the Fe-As interaction along the 
z-direction. Accordingly, the e-p coupling between this 
orbital and the A\ g mode would be expected to increase. 
This needs to be confirmed the future. 



to be explained within the conventional phonon mecha- 
nism. Nevertheless, the discovered dominance of the A\ g 
phonon mode in the e-p coupling may provide a clue to 
explain the large change of T c when the rare-earth el- 
ement is substituted by another one^ A possible role 
played by this type of phonons is the magnetoelastic ef- 
fect, since it is found that the magnetic coupling depends 
strongly on the z coordinate of the As atom^ It needs 
to be pointed out that after we have finished our calcula- 
tions we were informed by Egami that in BaFe2As2 the 
most important phonon mode is indeed of A\ g type ex- 
hibiting significant softening as a result of inelastic X-ray 
scattering measurements £0.. 

The dependence of e-p coupling on electron or hole 
doping needs to be investigated in future work, because 
the relative positions of the E g , B\ g flat bands with re- 
spect to the Fermi level change with doping. Therefore, 
the structure and strength of e-p coupling may change. 
Though these states dominate at the Fermi level as re- 
vealed in this work and other calculations, they have 
negligible e-p coupling in the present calculations. We 
have already shown that the Coulomb interaction is a 



IV. CONCLUSION 

Based on first-principles calculations, we have shown 
the similarity between the electronic structures of 
NaFeAs and other Fe pnictides. By using a tight-binding 
model, a global double degeneracy of bands along X-M 
and R-A directions has been revealed caused by the com- 
bination of site symmetry of Fe, As and the line symme- 
try of X-M in the BZ. The studies on the topology, met- 
ric property of FS and the generalized susceptibility x(l) 
demonstrate the paramagnetic state of NaFeAs to be un- 
stable against the Qj M -=(l/2,l/2,0) SDW state. Thus a 
SDW ground state with a possible structural distortion 
accompanied with symmetry breaking are predicted, a 
fact confirmed by recent experiments^ The possibility of 
a Qm related CDW can be completely excluded based on 
a resticted HFA criterion. A A\ g phonon involving only 
the z-direction vibration of As and Na has been detected. 
A possible enhancement of this specific e-p coupling due 
to the correlation effect has also been discussed. 
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